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ABSTRACT 

Numerical simulations predict that a large number of baryons reside in inter- 
galactic space at temperatures between 10^ — 10'' K. Highly-ionized metals, such 
as O VII and O VIII, are good tracers of this warm-hot intergalactic medium" , or 
WHIM. For coUisionally-ionized gas, the ionization fraction of each ion peaks at 
some particular temperature ("peak temperatures"), so different ions can there- 
fore trace the IGM at different temperatures. We performed a hydrodynamic 
simulation to study the metal distributions in the IGM. We then draw random 
lines-of-sight across the simulated region and synthesize resonance absorption line 
spectra in a similar way to simulating the Lya forest. By studying the distribu- 
tion functions of H- and He-like O, Si and Fe in a coUisionally-ionized IGM and 
comparing with semi-analytic results based on the Press-Schechter formalism, 
we find: (1) ions with higher peak temperatures (for instance, Fe XXVI) tend to 
concentrate around virialized halos, which can be well described by the Press- 
Schechter distribution, ions with lower peak temperatures are found both in small 
halos (such as groups of galaxies) and in filaments; (2) lower peak temperature 
ions are more abundant and should be easily observed; (3) peculiar velocities 
contribute a significant part to the broadening of the resonant absorption lines. 

Subject headings: intergalactic medium — large-scale structure of universe — 
quasars: absorption lines — X-rays: general 



1 Hubble Fellow 



- 2 - 



1. Introduction 

Cosmological N-body simulation plays an important role in studying cosmic structure 
formation and evolution. In this methodology, initial conditions are first set at some early 
epoch based on the history of the early universe and the linear growth theory of initial 
perturbations. N-body simulations are then used when cosmic structure grows into the 
non-hnear regime (for a review see Bertschinger 1998). The earhest cosmological N-body 
simulations involved only a simple system with a dissipationless dark matter component, and 
the evolution of such a system depends solely on gravity (see, e.g., Peebles 1970; Gott, Turner 
& Aarseth 1979). An important step was to incorporate the gas dynamics for baryonic matter 
(see, e.g.. Gen & Ostriker 1994; Kang et al. 1994. The evolution of this baryonic matter 
depends not only on gravity but also on heating and cooling processes. This step turned out 
to be quite successful; for example it reproduced the observed statistical properties of the 
Lya forest in spectra of high redshift quasars (see, e.g.. Gen et al. 1994; Zhang, Anninos & 
Norman 1995; Zhang et al. 1997; Machacek et al. 2000). In addition, numerical simulations 
based on the GDM models successfully fit the Press- Schechter distribution of galaxy clusters. 
Gombined with observations (see, e.g., Henry &; Arnaud 1991 for Einstein data. Burke et 
al. 1997 for ROSAT data and Henry (1997) for ASCA data), this method can be used to 
constrain cosmological parameters, such as fl and Ug (see, e.g.. Eke, Gole & Frenk 1996; Fan, 
BahcaU & Gen 1997). 

It is of particular interest to study the evolution of baryonic matter because it addresses 
an important question: why is there an apparent decrease in baryon density from high to 
low redshifts (see, e.g., Fukugita, Hogan & Peebles 1998)? Numerical simulations indicated 
that the absence of baryons at low redshifts represents a yet undetected IGM; the "warm-hot 
intergalactic medium" , or WHIM (Dave et al. 2000). For example. Gen & Ostriker (1999) and 
Dave et al. (2000) found that the average temperature of baryons is an increasing function 
of time, with most of the baryons at the present time having a temperature in the range 
of 10^ — 10^ K. Although invisible in the optical band, the highly ionized metals in WHIM 
gas might be detectable in UV and X-ray bands. Recent discovery of O VI absorption fines 
in the UV spectra of QSOs (see, e.g., Tripp, Lu & Savage 1997; Tripp, Savage & Jenkins 
2000) provided strong support for this, and implied a '^significant baryon reservoir'' at low 
redshifts. 

Since highly-ionized metals, such as O VII and O VIII, are also good tracers of WHIM 
gas at 10^ — 10^ K, in this paper we study the metal distributions in the IGM through 
hydrodynamic simulation. The idea is similar to that discussed in Perna & Loeb (1998), 
Hellsten, Gnedin & Miralda-Escude (1998) and Fang & Ganizares (2000); highly ionized 
metals would introduce absorption features in the X-ray spectrum of a distant quasar — 
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the X-ray Foresf . This paper has the following goals: (1) compare predictions of the X- 
ray forest from semi-analytic models to those from numerical simulations, this is important 
because each technique has its own strengths and weakness and by contrasting the two, we 
can determine which results arc robust; (2) compute predictions for a wider range of ions, 
including O VII, O VIII, Si XIII, Si XIV, Fe XXV, Fc XXVI; (3) explore the nature of the 
gas that gives rise to the absorption, including its thermal state and its location in clusters, 
groups or filaments. The method we use is in analogy to the method in simulating the 
Lya forest (see Zhang et al. 1997 for a complete discussion on formahsm), and is briefly 
outlined here: first, random lines-of-sight (LOS) are drawn across the simulated region, and 
temperature, density and velocity are obtained at each point along every LOS to calculate 
the optical depth due to the absorbing ions; every possible line feature is then examined 
with a certain opacity threshold and for each identified line, we fit a Gaussian profile to 
obtain its equivalent width (W), the Dopplcr 6-parameter and column density (N^). Here 
A^* = N{X^) is the column density of ion X*. Wc then obtain the X-ray forest distribution 
function (XFDF) by averaging a total of 5, 000 LOS and compare it with those obtained by 
the semi-analytic calculations in Fang & Canizares (2000). Similar work was also done before 
by Hellsten, Gnedin & Miralda-Escude (1998), based on a different ionization mechanism. 

This paper is arranged as follows: section 2 is a brief introduction to the numerical 
simulation, in section 3 we generate a simulated spectrum from the numerical simulation 
data and fit each individual line. We get the statistical properties of the fine distribution 
and compare them with the Press- Schechter formalism in section 4. Section 5 is a summary 
and discussion. 



2. Numerical Simulation 
2.1. Method 

We numerically simulated a cubic region of 100/i~^ Mpc, or 150 Mpc with h = 0.67 
on one side at z ~ 0. We model the dark matter with particles (as in a standard N-body 
code) but use the technique of adaptive mesh refinement (AMR) to follow the baryons. This 
method is designed to provide adaptive resolution combined with a modern, shock-capturing 
Eulerian hydrodynamics scheme. In this technique, a uniform mesh {^^parenf grid) covers 
the entire computational volume, and in regions that require higher resolution, a finer subgrid 
("c/izW grid) is added. By repeating this process recursively, high resolution can be achieved 
in small, high-density regions. To communicate between parent and child grids, boundary 
conditions go from coarse to fine, and the improved solution on the child grids is used to 
update the coarse parent grid. The advancing of the grid hierarchy is done automatically and 



-4- 



dynamically, so interesting features can be followed at high resolution without interruption. 
We refer to Bryan (1996) and Norman & Bryan (1999) for further details. 

In this simulation the dark matter particle mass was 5 x 10^^ Mq and the smallest cell 
size (in the high density regions) was A9h~^ kpc. To generate spectra, we created a uniform 
grid which contained 512 x 512 x 512 cells, with each cell having a size of ~ 195h~^ Kpc, 
or 19.5 km in velocity space. The simulation itself is a fiat low-density cold dark matter 
model (LCDM) with Qq = 0.3, = 0.7 and = 0.04. The amplitude of the initial power 
spectrum was fixed so that the rms overdensity in a spherical tophat with radius 8h^^ Mpc, 
as, is 0.9. We use a scale-invariant power spectrum with a primordial index of n = 1. The 
initial cosmological redshift is z = 30. 

2.2. Gas Distribution 

Figure 1 displays the projected baryonic density {Mproj, in unit of MqM^c"^) of the 
entire simulated region. Four panels are displayed with projected masses of Mproj > 10^°, 
10^^ 10^2 and 10^3 MqMpc-'^, respectively We can immediately distinguish four apparent 
features in Figure la: the red spots {Mproj > 10^"^ MqMpc^"^), which indicate regions with the 
highest densities; the yellow extended structures {Mproj ~ lO-*^^— 10^^ MqMpc^'^) surrounding 
the red spots; the green filamentary structures {Mproj ~ 10''^° — 10^^ MqMpc~'^) connecting 
the high density regions and the dark region {Mproj < 10^° M^Mpc"'^, for instance, on 
the top left corner) indicating regions with the lowest densities. By looking from panel (a) 
through (b) to (c) and (d) we find that as the cut-off Mproj increases, most of filaments 
begin to disappear and high density regions start to concentrate only on several isolated red 
spots (for example, see panel c and d). These highly mass-concentrated regions indicate that 
collapsed, virialized objects, such as groups and clusters of galaxies, have formed there. Since 
overlap of these massive systems is rare along the projected direction, a projected surface 
density of 10^^ or 10^^ MqMpc'"^ roughly corresponds to a group or cluster of galaxies with 
M ~ 10^^ or 10^^ Mq, respectively. 

Figure 2 shows the emission-weighted temperature distribution of the simulated region. 
The four plots display the temperature distributions in four ranges: T > 10'', 10^, 10^ and 
10^ K. The temperatures below about 2 x 10^ K are not real because the simulation does 
not include an ionizing radition field which would heat the gas up. The highest temperature 
in our simulation is 8.09 x 10^ K. As shown by figure 2a, 2b, 2c, most of the simulated 
regions are filaments with 10"^ K < T < 10^ K. The intersections between filaments show 
higher temperatures with T > 10^ K. Comparing Figure 2 with Figure 1, we find that high 
temperature regions correspond to high density regions. For instance, comparing panel c in 
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both figures, a projected surface density of 10^^ MqMpc ^ implies a group of galaxies with 
M ~ 10^2 Mq and T ~ 10^ K, as indicated in Figure 2c. 

To ease the load of numerical computing, no atomic and radiative processes are taken 
into account here, so high gas temperature is mainly due to adiabatic compression and shock 
heating. Gas heating and cooling processes will play an important role in the cold gas (and 
stars) in galaxies, but should play only a minor role for hot gas. Dave et al. (2000) compared 
an AMR simulation including only adiabatic process (model Bl in their paper) with an AMR 
simulation including several gas heating and cooling processes (model B2). These processes 
include H, He and metal cooling, photoionization and star formation & feedback. They 
checked the baryon fraction of WHIM gas and found that both numerical simulations gave 
consistent results at 2; = {^wHiM/^b = 0.32 for model Bl, and 0.3 for model B2). By 
comparing the two AMR simulations with several other numerical simulations they found 
that gravitational shock heating of gas faUing onto large scale structures is the dominant 
heating mechanism for WHIM gas. Most of the higher temperature IGM (T > 10^ K) is 
associated with hot intracluster gas (see Figure 1 and 2). The major cooling mechanisms for 
intracluster gas are thermal bremsstrahlung radiation and X-ray line emission. The cooling 
time scale is approximately given by (Sarazin 1988) 



which is longer than the Hubble time. Here Hp is the proton density and typically in the 
order of 10~^ cm~^. Based on the above considerations, we find that excluding gas heating 
and cooling in our simulation is a reasonable approximation for an IGM with T > 10^ K. 



Little is known about the metallicity of the IGM. The metal abundance could be as low 
as IO^^Zq in the low column density Lya forest systems (Lu et al. 1998) or as high as IZq 
in the interstellar medium (ISM) of some early-type. X-ray luminous galaxies (Matsushita, 
Ohashi, & Makishima 2000). Since we will compare numerical simulation with the semi- 
analytic model in our previous paper (Fang & Canizares 2000), we assume uniform metalicity 
and use the metallicities observed from the intracluster medium, e.g., 0.5Zq for oxygen, 
sihcon and O.SZq for iron (Mushotzky et al. 1996). Table 1 lists parameters for the ions and 
lines considered here. 

To investigate X-ray absorption in the spectrum of a background source, we need to 
understand the distribution of ionization states of metals in the IGM. To obtain the ion 
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density, we assume the gas is in coUisional equilibrium and we adopt the ionization fractions 
from Mazzotta et al. (1998). We select helium- and hydrogen-like oxygen, silicon and iron. 
The peak temperatures of these ion species range from 10^ K to 10^ K (see Figure 3), which 
probe different temperature regions of the IGM. 

Figure 4 through Figure 9 display the projected column density, distributions of all 
the six ions. In each figure, the ion column density distributions are plotted in four different 
panels with A^* > 10^^, 10^^, 10^^ and 10^^ cm^^, respectively. The color bar on top of each 
panel indicates the logarithm scale from the minimum to the maximum column densities in 
each panel. 

By examining the six figures, wc find that the distributions of these ions are closely 
related to the density (Figure 1) and tempcraturc(Figure 2) distributions. In general, column 
density distributions of metals follow the projected density distribution in Figure 1: regions 
with higher baryon densities also have higher metal densities. Due to the fact that different 
ions have different peak temperatures, there are differences among the distributions of these 
ion species. For instance. Figure 2d shows that many of the filaments that connect high 
density regions are shock-heated to temperatures above 10^ K. Since O VII has a peak 
temperature between 3 x 10^ K and 2 x 10^ (Figure 3), we find that a large amount of O VII 
exists in filaments with a column density even higher than 10^^ cm~^ (Figure 4c). On the 
other hand, Fe XXVI has a peak temperature of over 10*^ K. In Figure 9 wc find that almost 
no Fe XXVI exists in filaments, and it all concentrates around those collapsed intersections 
that appear in Figure 1. By comparing these two extremes, we find that generally Fe XXVI 
follows the distribution of virialized objects while most of O VII exists in filaments, and the 
distributions of other ion species vary between those of Fe XXVI and O VII. 



3. Spectral Analysis 

The methodology we adopt is straightforward: a random linc-of-sight (LOS) is drawn 
across the simulated region; temperature, baryon density and velocity are obtained at each 
point along the LOS; the ion density is calculated based on the metallicity and ionization 
fraction; the synthesized spectrum is then obtained from the optical depth along the LOS; 
finally, by fitting the resonant absorption lines in the "simulated' spectrum we obtain ion 
column densities {N'^) and the Doppler parameters (b). This method is similar to the method 
used in the Lya forest; we refer to Zhang et al. (1997) for a complete discussion and formal- 
ism. 
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3.1. Gas Distribution along the LOS 

To investigate the gas distribution along the LOS, we draw random hnes across the 
simulated region. Each LOS has 512 points. The separation between each point is 0.195/i^^ 
Mpc, or IS.Z = 6.5 x 10^^ in redshift space. This roughly corresponds to an energy resolution 
of 0.04 eV around 0.6 keV (oxygen line region), or 0.44 eV around 6.7 keV (iron line 
region). We obtain temperature, peculiar velocity (coherent or turbulent) and gas density 
at each point. For demonstration, we show their distributions along one LOS in Figure 10: 
Figure 10a shows the temperature distribution. Figure 10b gives the velocity distribution 
and Figure 10c gives the baryon number density distribution. 

The dotted line in Figure 10c is the mean baryon density of the universe 

(p) = 1.124 X IQ-^Qbh^ cm-^ = 2.136 x 10"^ cm"^ (2) 

Here 0,b is the baryon density and we adopt a value of Qbh^ — 0.019 from the standard theory 
of big bang nucleosynthesis (BBN) (see, e.g.. Buries, NoUett & Turner 2000). Most of the 
regions along this random LOS are underdense regions, with several peaks indicating over- 
dense regions. The largest overdense region occurs at a distance of approximately 50 Mpc, 
with 6p = p/{p) ~ 20. Given the temperature T ~ 10^ — 10^ K (see the corresponding 
position in Figure 10a), this region might correspond to a small cluster or group of galaxies. 

We notice that one region (at ~ 110 Mpc) shows both a very high temperature (2— 5x 10^ 
K) and a high peculiar velocity, while this region is still underdense. This indicates that the 
LOS passed near a collapsed structure. The shock waves from collapsed structures can 
propagate deeply into low density regions, which cause a strong divergence of the peculiar 
velocity, and in turn a high temperature in this underdense region. 

The density peaks of different ion species are closely related to peaks of their ionization 
fraction curves (see Figure 11). For instance, the peak temperature of Si XIV, Fe XXV and 
Fe XXVI are all above 10^ K, and since the only region with such a high temperature along 
this LOS is at ~ 110 Mpc (see Figure 10a), these ion species have just one density peak 
around that region (the bottom three panels of Figure 11), although this region is below the 
mean density of the universe. 

3.2. Simulated Spectra 

Considering the X-ray spectrum of a background quasar located at redshift Zo, the 
frequency of a photon emitted at will he u = uq^I + zo)~^ when the photon reaches the 
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observer. The optical depth at frequency v (observer frame) is (Spitzer 1978) 

r^" (11 
r(i/) = j m (z) ai (i/') —dz (3) 

where ni{z) is the comoving number density of absorbing ion i at a rcdshift of z and ai is 
the ion absorption cross section at frequency f'. Besides the Hubble velocity, the peculiar 
velocity (both coherent and turbulent) also contributes to the line broadening and produces 
a shift at the hne center. Assuming a pecuhar velocity of Vpec along the line-of-sight (LOS), 
the frequency at redshift z is 

jy'^u{l + z){l+^-^) (4) 

To the first order (vpec/c) approximation the optical depth at frequency u (observer frame) 
is (Zhang et al. (1997)) 




(5) 



Here = {ne'^ /nipc)/ is the absorption cross section at the line center (/ is the upward 
oscillator strength) and 6r is the velocity dispersion due to the thermal motion. 

By integrating over a distance of 150Mpc (^o ~ 0.034 in our cosmological model) we 
synthesize the absorption spectrum. The solid lines in Figure 12 show the transmission, 
F = exp (— r), from O VII to Fe XXVI (top to bottom), where r is calculated from equation 
(5). The translation from redshift {z) to wavelength (A) is simply 

(6) 

where Aj is the ion's rest wavelength. 

Comparing Figure 11 and Figure 12 we can see that absorption features are created 
by density peaks. Along this random LOS we find that most of the lines are very weak 
(r < 0.01), except for some O VII lines, where the largest optical depth at the line center 
is almost t ^ 2. The densities of Fc XXVI are so small (less than 10~^^ cm~^) that no 
absorption feature is displayed in the bottom panel of Figure 12. 



3.3. Line Fitting 



To exactly measure absorption lines we need to fit the synthesized spectrum with a 
Voigt line profile. Numerically such a procedure is difficult to achieve. However, the line 
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profile can be well represented by a Gaussian line profile in our situation (see Zhang et al. 
1997 for discussion). To simplify we adopt the threshold-has ed'' method of spectral analysis 
developed by Miralda-Escude et al. (1996) and Hcrnquist et al. (1996). In this method, we 
follow the flux of a simulated spectrum. We flrst define a threshold Ft. An absorption feature 
is identified whenever the fiux drops below (1 — Ft) at wavelength Ai and rises above 1 — Ft 
at wavelength A2 where A2 > Ai. In case of line blending between Ai and A2, different lines 
are identified whenever the slope of transmission curve changes sign and fitting is performed 
for each identified line. Here we see that Ft actually represents the minimum central optical 
depth of lines which can be resolved by this method. The line equivalent width (W) is 
calculated by integrating (1 — F^) from Ai to A2. Taking into account redshift, the equivalent 
width is given by 

The Doppler 6-parameter is defined as the 6-parameter of a Voigt profile which gives 
the same equivalent width with internal width A A = A2 — Ai between Ai to A2. Following 
Spitzer (1978) we have 

■ (8) 



2AoF(ro) 

Here tq is the optical depth at the line center, and 



F(to)-^ [1 - exp (-Toe-^')] (9) 



is the curve-of-growth function. 

The ion column density is given by 

= ^ (10) 
c^oAo 

Here ctq is the absorption cross section at the line center. 

We need to understand the effect of threshold on line identification before we analyze 
spectral properties. Equation (10) shows that for a given threshold, there exists a minimum 
resolvable column density: the ^' threshold-has ed'' method cannot resolve a line with flux above 
this threshold. The lower the threshold is, the lower the minimum column density is. To 
understand this effect, we fit the synthesized O VII and O VIII spectra for the 5000 random 
line-of-sight (LOS), with different thresholds ranging from Ft — 10~^ to Ft — 0.8. In Table 2, 
we list thresholds (the first column) and the minimum resolved column densities for O VII 
(the second column) and O VIII (the third column). As we discussed in Fang & Canizares 
(2000), the minimum detectable column density for oxygen, even with the future mission 
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Constellation-X, is ~ 10^^ cm""^, which means that a high threshold of Ft — 0.5 — 0.8 is more 
reahstic from an observational point of view. However, in order to understand the low density 
regions with low metal column densities and make our sample theoretically complete, in the 
following discussion we select the lowest threshold. Ft = 10"^, to study spectral properties 
statistically. This also helps in comparing the numerical simulation with the semi-analytic 
calculations we performed in Fang & Canizares (2000). 

Figure 12 shows six spectra from the sample LOS in section §3.2 ranging from O VII to 
Fe XXVI. The solid lines are the spectra from the numerical simulation and the dashed lines 
are from the fitted data by the ^'threshold-based^ method. Here Ft = 10^^ (the dotted line). 
We identify five O VII, four O VIII, five Si XIII, two Si XIV, one Fc XXIV and zero Fe XXVI 
absorption lines. 

The line parameters are listed in Table 3. The first column is the ion species; the second 
column is the absorption line numbers; the third column gives the line redshift; the fourth 
column is the column density; the fifth column is the equivalent width; the final column 
gives the Doppler 6-parameter. We can see that all fines are identified, even including those 
which are blended. For example, in the transmission spectrum of Si XIII in Figure 12, lines 
1 and 2 are blended with each other. Since the sign of the curve changes from line 1 (plus) 
to line 2 (minus), two lines are identified (sec Table 3) with a separation of Az = 3.3 x 10~^. 
Table 3 shows that the equivalent widths of most absorption lines are so weak (W < O.leV) 
that they are practically undetectable. However, this table does show some O VII lines 
with W > O.leV, which might be detectable by future missions, such as Constellation- X. 
The Fe XXV line shows a large b value, indicating that this line is generated in a high 
temperature region. 

4. Spectral Properties 

In this section we discuss the statistical properties of the absorption lines identified 
in our numerical simulation. We draw a total of 5000 random LOS across the simulated 
region. For each LOS we generate a simulated spectrum and fit it by the ^'threshold-based^ 
method. Table 4 lists the maximum and minimum value of column density and the Doppler 
6-parameter for each ion for Ft — 10~^. In the following sections we will discuss column 
density and the Doppler ^-parameter distributions of these lines. 
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4.1. Column Density Distribution 

In analogy to the Lyman-o; forest system, we define d'^P/dN^dz as the number of 
absorption systems along the line of sight with a column density between and + dN'^ 
per unit redshift. This is the so-called "X-ray forest" distribution function (XFDF). The 
calculation of XFDF is straightforward: for each ion, the column density is divided into 
several bins between 10^^ and 10^^ cm~^; absorption lines are counted in each column density 
interval and then divided by Az. Here Az = 3.4 x 10~^ is the size of the simulated region 
in redshift space. The lower limit of the column density is selected arbitrarily; although a 
lower limit of 10^^ cm~^ is undetectable, we still include it to make our sample theoretically 
complete; the upper limit is selected by noticing that no ions shows a column density larger 
than 10^^ cm~^ in Table 4. 

Figure 13 displays the distribution functions from O VII to Fe XXVI. The solid lines show 
all the helium-like ion species, while the dashed lines show all the hydrogen-like ion species. 
The main feature seen in this figure is that all curves display a sharp cutoff somewhere 
between 10^^ and 10^'^ cm^^. The reason for this apparent drop at high column density end 
is two fold: first, systems with high column density generally correspond to high temperature 
regions (for example, compare Figure 1 with Figures 4 through 9), and at high temperature 
the ionization fraction drops quickly and so does the ion column density; second, as we 
emphasized in section §2.3, high column density systems are closely related to virialized 
objects, such as groups and clusters of galaxies, and the Press- Schechter theory predicts 
that the spatial density of virialized objects decreases exponentially as their masses (and 
temperatures) increase. Based on this explanation, the maximum ion column densities (iV^^^ 
in Table 4) should come from the centers of virialized objects which have temperatures close 
to the peak temperature of that ion. For instance, a rough estimation shows that for a galaxy 
cluster at temperature 5 x 10^ K, the central column density of Fe XXV is ~ 4.7 x 10^^ cm~^, 
which is about the same order as the value shown in Table 4. 

Figure 13 also shows that the overall amplitude of the XFDF is a decreasing function 
of the ions' peak temperature: ions with low peak temperatures display higher amplitudes 
which implies more absorption lines (see table 4) . This is not only due to the differences of 
the metal abundances, but also due to the fact that most of the baryons have temperatures 
ranging from 10^ — 10^ K and the baryon fraction decreases as the temperature increases 
(Dave et al. 2000). This imphes that the low peak temperature ions such as O VII and 
O VIII are more abundant and should be easily observed. 
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4.2. Press-Schechter Distribution 

Returning to Figure 1 and 2 we find that in general mass and temperature distributions 
follow each other: the low-density regions, such as filaments, show temperatures between 
10^ and 10^ K; the high-mass regions, such as intersections between filaments, show higher 
temperatures with T > 10^ — 10^ K. High temperature and high density indicate that 
collapsed, virialized objects, such as groups and clusters of galaxies, might be formed in 
these intersection regions. As we emphasized before, due to the different peak temperatures 
of the ion species, these ions actually probe the different temperature regions of the IGM. 
For instance, O VII and O VIII, which have peak temperatures around 10® K, should follow 
the distributions of those filamentary structures; on the other hand, Fe XXV and Fe XXVI, 
which have peak temperatures over 10^ K, should follow the distributions of those collapsed, 
virialized objects. Figures 4 through 9 confirm this point. 

Based on the formalism in Fang & Canizares (2000), we compare our numerical sim- 
ulation with the analytic model. The Press-Schechter function (Press & Schechter 1974) 
predicts the spatial distribution of virialized objects. We would expect that the analytic 
model and the numerical simulation should give consistent results for the distribution of 
the high-peak-temperature ions, such as Fe XXV and Fe XXVI, because these ions follow 
the distribution of virialized objects; on the other hand, we would not expect the low-peak- 
temperature ions, such as O VII and O VIII, would have the same behavior because most of 
these ions are distributed in filaments. 

Before we compare the ion column density distributions, we first need to compare the 
spatial number density distribution of the virialized objects identified in our numerical simu- 
lation with that predicted by the Press-Schechter function. To identify the viriafized objects 
in the simulated regions, the spherical overdensity algorithm (Cole & Lacey 1996) is applied. 
In this algorithm, peaks in the density distribution are identified as the center of virialized 
objects; around each point, a sphere grows until the mean density inside the sphere reaches 
some critical value, A^. The mass center of this sphere is computed and used to grow a new 
sphere. This procedure is iterated to convergence. If the center-to-center distance of two 
identified objects is less than 3/4 of the sum of their virial radii, we define these two objects 
as overlapping and the less massive object is removed (see Bryan & Norman 1998). 

Figure 14 displays the comparison between the numerical simulation (the filled square) 
and the analytic model (the sold line) for the number density of virialized objects. The 
analytic calculation is based on the LCDM cosmological model we adopt in the numerical 
simulation. The agreement is generally good, although the low-mass objects identified by 
the numerical simulation are systematically fewer than predicted. This is due in part to the 
fact that those low-mass objects are poorly resolved. However, there is also evidence in many 
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numerical simulations indicating that the disagreement seems to be a general shortcoming 
of the PS method (sec Laccy & Cole 1994; Gross et al. 1998). Governato et al. (1998) found 
that at o"8 = 0.7 or higher, the PS approach overestimates the number of small halos with 
M < lO-'^^/i^^ Mr.-j, which is consistent with our result. This deficit has been associated 
with the merger events not accounted for by the PS formalism (Cavaliere & Menci 1997; 
Monaco 1997). By studying the clustering of less massive halos Sheth & Tormen (1999) 
revised the PS function with a bias factor and obtained a good fit to numerical simulations 
for small halos as well as large halos. In our case the agreement in the high-mass region 
(M > 10^^/i~^ Mq) is very good. 

The techniques used to calculate the XFDF have been described in Fang & Canizares 
(2000), so we only briefly summarize them here. The mass of a viriahzed halo is 



Here pcru is the critical density of the universe, r^r is defined as the virial radius of a 
spherical volume within which the mean density is AcPcrit, and we adopt a conventional 
value of Ac = 200. The mass-temperature relationship from Eke, Cole & Frenk (1996) is 
used to convert between temperature and mass. Finally, to obtain the ion column density 
distribution within the halo we use the /3-model of intracluster gas density distribution 
(Sarazin 1988) and assume a gas fraction of fgas — 0.2. 

Since ionization fractions are very sensitive to temperature, we also investigate the 
assumption that virialized halos are isothermal. In Figure 15 we plot the normalized radial 
temperature profile for several identified clusters in the simulation. We find the overall radial 
temperature profile can be well represented by (the solid line) 



Here Tyij. is the virial temperature of clusters. The computed ion column density using the 
PS formalism takes account of this temperature variation along the integral path. 

Figure 16 shows the XFDF predicted from the Press- Schechter approach for all the 
six ions (solid line in each plot). The filled-dotted line displays the distributions from the 
numerical simulation. The most striking feature is that the predicted distributions fit the 
simulation very well at high column density regions, specifically at > 10^^ cm^^ for 
both oxygen species and at Ni > 10^^ cm~'^ for sihcon species and Fe XXV. For Fe XXVI, 
both distributions agree very well in the whole column density range (10^^ — 10^^ cm~^). 
This result can be well understood in the framework we describe at the beginning of this 
section: high-peak-temperature ions, such as Fe XXVI, follow the distribution of collapsed. 



(11) 




(12) 
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virialized halos well described by the PS formalism. For low-peak-temperature ions, the high 
column densities also come from virialized structures, while those with low column density 
are distributed in those filaments which can not be described by the PS distribution. Thus 
except for Fe XXVI, at low column density the PS distributions are systematically lower 
than in the numerical simulation. 

For high column density O VII and O VIII (iV > 1 - 2 x 10^^ cm-^) we find that the PS 
distributions are slightly higher than the distributions from the numerical simulation. This 
is because the high column density oxygen lies in virialized halos with a temperature around 
10^ K, which corresponds to a mass between lO^^ ^ IQI^Mq. As we find in Figure 14, the 
PS formalism overestimates the number of small halos, so the XFDFs for O VII and O VIII 
at > lO^^cm"^ are also overestimated. 



Dopplcr 6-paramctcr governs the line width and can be affected by several factors, 
such as gas temperature (thermal motion), the Hubble flow and the velocity structure (see 
equation [5]). We follow the approach used in the Lya forest studies (Bryan, Machacek, 
Anninos, & Norman 1999; Bryan & Machacek 2000) and distinguish two sources of fine 
broadening: bx, the thermal broadening; bprof, a convolution of the density profile with the 
peculiar velocity. The total 6-parameter is 



Figure 17 shows the 6-parametcr distribution function for six ions (the filled circles). All 
six ions exhibit a low cutoff to the Dopplcr parameter. This cutoff is also listed in Table 4. 
The cutoffs arc primarily due to the cutoff temperatures for each ionization fraction (see 
figure 3). For instance, the minimum temperature we used for O VIII is T^ut = 6.31 x 10^ K, 
so 



which is about the same value listed in Table 4. However, this also imphes the fact that the 
thermal motion may play a dominant role in line broadening {bprof ~ and b ^ br)- 

To understand the role of different line broadening mechanisms, we fit the 6-parameter 
distribution with a Gaussian (the dashed fine in Figure 17) 



4.3. Doppler Pcirameter Distribution 




(13) 




(14) 



/ocexp [-(6 -6)2/2(72] 
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Here b is the mean velocity and a is the standard deviation of a Gaussian curve. Both 
parameters are hsted in the column 2 & 3 of Table 5. By comparing b for different ions we 
can find that 



1 



6oc(^^J (16) 

Here Tpeak is the peak temperature for different ions shown in Figure 3. This confirms 
that most absorbing ions come from a very narrow temperature range around their peak 
temperature Tp^ak- 

However, we also find that there exists a non-Gaussian tail in each plot. This tail can 
be well fitted by a lognormal distribution (the solid line) 

/ocexp[-log^(6/&^n)/2af„] (17) 

Here 6/„ and o"/„ are the mean velocity and the standard deviation of this lognormal dis- 
tribution. Both fitting parameters are listed in columns 4 and 5 of Table 5. We find this 
distribution naturally cuts off at low 6, and also provide an excellent fit to the long tail. This 
non-Gaussian tail implies that the Hubble flow and peculiar velocity contribute a signiflcant 
part to those ions with high b values. 

To associate the Doppler 6-parameter with column density for each line, we show the 
scatter plots of b vs. N'^ for all six ions (Figure 18). Each plot displays a left border and 
a lower cutoff. The left border is an artifact of selecting a threshold Ft in synthesizing 
spectrum. For instance, in the case of O VII, given Ft — 10~^, equation [10] imposes (the 
dotted line in Figure 18) 

IQi 



6< 2250 (— ^) kms-^ (18) 



To understand the lower cutoff b^i^, we need to understand where the lines with the 
lower cutoff come from. From Figure (18) we know that for lines with a given column density 
iV*. they have different Doppler-6 values. These lines come from either virialized halos (the 
intersection regions) or filaments. Since filaments are shock-heated (see, e.g., Dave et al. 
2000) and have a much larger bprof, we assume that the lines with the lowest Doppler-6 

1 /2 

values (bmin) come from virialized halos, and for these lines bmm ^ bx oc Tj^. where T^ir is 
the virial temperature of the halo. However, even these lines could come from different halos 
with different virial temperatures: they could come from either the outskirts of a larger halo, 
or the central region of a smaller halo. According to the mass-temperature relationship we 
know that a smaller halo has a lower virial temperature and so has a smaller b value. We 
conclude that those lines with lower cutoffs (bmin) come from the center of virialized halos. 
In Figure 18 we plot br vs. the central column densities of all virialized halos (the solid 
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lines), based on our analytic model. We find that this gives a quite good fit to the lower 
cutoff in all six plots. The turnover at higher b is due to the fact that the ionization fractions 
fall at higher temperatures and so do the central ion column densities. 

5. Summary and Discussion 

In this paper we have performed a hydrodynamic simulation based on a LCDM model. 
By drawing random LOS and synthesizing the X-ray spectrum for each LOS, we obtained 
the XFDFs for hydrogen and hclium-like O, Si, Fe. We have investigated the column density 
distributions of highly- ionized metals by comparing the numerical simulation to the analytic 
method described in Fang & Canizares (2000). We also discuss the distribution of Doppler 
6-parameters. Our main conclusions can be summarized as follows: 

1. Ions with low peak temperatures such as O VII and O VIII trace small halos (galaxy 
groups) and filaments, while high peak ions such as Fe XXV and Fe XXVI are only 
observed in the largest objects (galaxy clusters). For all ions, the strongest signal come 
from virialized halos; however, filaments show up as lower column density absorbers 
(see Figure 4 to 9). 

2. The low peak temperature ions are more abundant than ions with high peak temper- 
atures. This can be understood from numerical simulations that most of the baryons 
have temperatures ranging from 10^ — 10'' K and the baryon fraction decreases as the 
temperature increases. This implies that from the point view of observation, ions with 
low peak temperatures such as O VII and O VIII should be more easily observed. 

3. We compared the column density distribution of absorbers predicted by the simulations 
against an analytic model developed in Fang & Canizares (2000) and found that 
the numerical simulation confirms our predictions: ions with high peak temperatures 
(T > lO^K) follow the distributions of collapsed halos, which can be described by the 
Press-Schechter distribution; for ions with low peak temperatures, those lines with high 
column density come from virialized structures, while those with low column density 
are distributed in diffuse filaments which are not well described by the PS distribution. 

4. Higher excitation ions produce lines which arc generally wider (as determined by the 
Doppler 6-parameter) than low excitation lines, although there is a wide distribution 
in all cases. This distribution is well fit by a log normal distribution. The sharp lower 
cutoff of this distribution is provided by the fact that there is a minimum temperature 
required to create the ion in the first place. Peculiar velocities in the gas also contribute 
significantly, providing a non-Gaussian tail at the high-6 end. 
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The largest uncertainty in spectral synthesis comes from the metallicity: we did not 
include any gas processes in our numerical simulation. Although these gas processes are not 
important in heating and cooling gas with T > 10^ K, the IGM metals enriched by these 
processes are crucial in UV/X-ray observations: it is these metals which produces detectable 
features. We assume a uniform metallicity in generating the synthesized spectrum, however, 
we would expect that metallicity should be determined by the metal enrichment history of 
the IGM, which in turn depends on a variety of sources: the UV/X-ray background radiation, 
supernova input, etc, and these are likely to be different in low and high density regions. 
It is important to incorporate all these processes to obtain a self-consistent result of metal 
distributions and ionization structures in the IGM. This is the goal of our future numerical 
simulations. 
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Table 1. Ion Parameters 



Ion 


Energy 


Wavelength 


Oscillation 


MetaUicity^ 


Solar 




(keV) 


(A) 


Strength 




Abundance'' 


VII 


0.57 


21.6019 


0.696 


0.5 


8.53 X 10"^ 


VIII 


0.65 


18.9602 


0.416 


0.5 


8.53 X 10"^ 


Si XIII 


1.87 


6.6480 


0.757 


0.5 


3.58 X 10-5 


Si XIV 


2.01 


6.1822 


0.416 


0.5 


3.58 X 10-5 


FeXXV 


6.70 


1.8504 


0.798 


0.3 


3.23 X 10-5 


Fe XXVI 


6.97 


1.7798 


0.416 


0.3 


3.23 X 10-5 



^Relative to solar abundance Zq. 
^Relative to hydrogen number density. 
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Table 2. Minimum Column Densities for O VII and O VIII 



Threshold (Ft) 


VII (cm-2) 


VIII (cm-2) 


10-^ 


1.43 X 10^° 


2.75 X 10^° 


10-1 


1.34 X 10^^ 


4.28 X 10^^ 


0.2 


2.35 X 10^^ 


8.37 X 10^^ 


0.5 


6.82 X lO^^ 


1.96 X 10^^ 


0.8 


1.13 X 10^^ 


3.75 X 10^5 
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Tablc 3. Resonant Absorption Lines 



Ion 


Line 


z 




N (cm-2) 


EW (el/) 


b {km/s) 


VII 


1 


8.57 X 10" 


-3 


6.40 X 10^^ 


2.72 X 10-^ 


76.05 







i.i / X iU 


-2 


i.U4 X iU 


0.4D X iU 


/ y.io 




3 


2.19 X 10" 


-2 


3.84 X 10^^ 


2.42 X 10-2 


23.00 




4 


2.46 X 10" 


-2 


1.74 X lO^-^ 


1.23 X 10-3 


44.93 




5 


2.94 X 10- 


-2 


7.76 X 10^2 


5.41 X 10-^ 


16.88 


VIII 


1 


1.13 X 10" 


-2 


1.23 X 10^^ 


5.42 X 10-3 


128.90 




2 


1.19 X 10" 


-2 


1.05 X 10^^ 


4.57 X 10-3 


53.98 




3 


2.46 X 10" 


-2 


1.03 X 10^3 


4.37 X 10-^ 


55.61 




4 


2.64 X 10" 


-2 


2.83 X 10^2 


1.20 X 10-^ 


89.50 


Si XIII 


1 


1.08 X 10" 


-2 


5.75 X 10^2 


4.63 X 10-^ 


56.14 




2 


1.11 X 10" 


-2 


6.36 X 10^2 


5.11 X 10-^ 


65.85 




3 


1.18 X 10" 


-2 


9.45 X 10^2 


7.57 X 10-^ 


42.65 




4 


2.46 X 10" 


-2 


1.13 X 10^2 


8.73 X 10-5 


60.77 




5 


2.64 X 10" 


-2 


2.27 X 10^2 


1.73 X 10-^ 


75.55 


Si XIV 


1 


2.47 X 10" 


-2 


1.41 X 10^2 


5.98 X 10-5 


267.95 




2 


2.62 X 10" 


-2 


6.22 X 10^1 


2.63 X 10-5 


90.47 


FeXXV 


1 


2.44 X 10" 


-2 


3.24 X 10^2 


1.93 X 10-6 


220.38 



^Prom a sample LOS, see Figure 12. 
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Table 4. Line Properties 



Ion 


Line Number 


' mm 


' max 














{km/s) 


{km/s) 


VII 


10222 


1.43 X 10^° 


8.39 X 10^^ 


13.43 


906.70 


VIII 


7078 


2.75 X 10^° 


4.57 X 10^^ 


25.73 


829.26 


Si XIII 


6246 


4.67 X 10^° 


1.30 X 10^6 


25.81 


697.71 


Si XIV 


2960 


1.34 X 10^^ 


1.13 X lO^*' 


39.30 


454.34 


Fe XXV 


1192 


2.86 X 10^1 


2.56 X 10^6 


56.36 


483.18 


Fe XXVI 


194 


8.51 X 10^^ 


5.51 X 10^^ 


68.63 


356.05 
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Table 5. Doppler Parameter Distributions 



Gaussian Log normal 



Ion b a bin crv. 



VII 


47.17 


23.06 


38.52 


0.23 


VIII 


54.84 


22.47 


51.97 


0.19 


Si XIII 


62.29 


30.48 


50.03 


0.22 


Si XIV 


94.18 


34.85 


84.02 


0.17 


Fe XXV 


96.14 


29.63 


89.51 


0.14 


Fe XXVI 


124.10 


30.76 


116.17 


0.11 
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id) 13.0 Hsa^i:-:'.::".: 1 S.7 




M > 10'^ M^pc 



Fig. 1. — The projected surface density {Mproj, in unit of MqMpc~'^) of the entire simulated 
region. Four panels are displayed with projected masses of Mproj > 10^°, 10^^, 10^^ and 
10^^ MqMpc~^, respectively. The color bar on top of each plot indicates the color scale in 
logarithm. 




T > 10' K T > 10" K 

Fig. 2. — The emission- weighted temperature distribution of the simulated region. The four 
plots display the temperature distributions in four ranges: T > 10*, 10^, 10^ and 10^ K. 
The highest temperature in our simulation is 8.09 x 10^ K. The temperatures below about 
2 X 10^ K are not real. 
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Temperature (K) 

Fig. 3. — Ionization fractions from O VII (solid curve) to Fe XXVI (long-dashed curve). 





Fig. 4. — The projected density (column density, iV*) distributions of O VII, plotted in four 
different panels with N'^ > 10^^, 10^^, 10^^ and lO"*^^ cm~^, respectively. The color bar on 
top of each panel indicates the logarithm scale from the minimum to the maximum column 
densities in each panel. 



(a) LS.O 



MJiSM 17.0 




(t) 13.0 



isri(o vnr) > lo'^ cm- 

(c) 14.0 ^■nstS'Vf.v. ?msa 17.0 





N,(0 VIII) > 10" cm- 

id) 15.0 Hsa^i:-:'::".: 1 7.0 
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Fig. 5. — The projected column density distributions of O VIII. 
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Fig. 6. — The projected column density distributions of Si XIII. 
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Fig. 7. — The projected column density distributions of Si XIV. 



(a) l^.D Hnsi^vf.vv. 16.4 (t) 13.0 Hsa^::v.::".: rsf^ 1 6.4 




Wi(Fe XXV) > 10'* cm"' N,{Fe XXV) > 10" cm^ 

Fig. 8. — The projected column density distributions of Fe XXV. 
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Wi(Fe XXVI) > 10'' cm^ N,{Fe XXVI) > 10" cm^ 

Fig. 9. — The projected column density distributions of Fe XXVL 
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Fig. 10. — (a) Temperature, (b) peculiar velocity (coherent or turbulent) and (c) gas density 
distributions along one LOS. The dotted line in (c) is the mean baryon density of the universe. 
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Fig. 11. — The number density distributions along one LOS for six ions. 
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Fig. 12. — Synthesized absorption spectra. The sohd hnes are the spectra from the numerical 
simulation and the dashed lines are from the fitted data by the ''Hhreshold-based' method. 
Here Ft — 10~^. Dotted lines indicate 1 — Ft. 
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Fig. 13. — The distribution functions from O VII to Fe XXVI for Ft = 10 ^. 



-39- 




1 ? 1 14 

10^^ ^qIO ^qI^ 

M (Me) 



Fig. 14. — The comparison between the numerical simulation (the filled square) and the 
analytic model (the sold line) for the number density of virialized objects. The analytic 
calculation is based on the LCDM cosmological model we adopt in the numerical simulation. 
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Fig. 15. — The normalized radial temperature profile for several identified clusters in the 
simulation. The solid line is a fit to simulated profiles, T^j^ and r^ir are the virial temperature 
and radius of clusters. 
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Fig. 16. — The XFDF predicted from the Press-Schechter approach (sohd hnes) and numer- 
ical simulation (fiUed-dotted lines). 
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Fig. 17. — The 6-parameter distribution function for six ions (the dotted hues), fitted with 
Gaussian distributions (the dashed hues) and log normal distributions (the sohd lines). See 
Table 5 for fitting parameters. 
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Fig. 18. — The scatter plots of b vs. A^* for all six ions. Each plot displays a left border, 
which gives a maximum b value (bmax), and a lower cutoff, which gives a minimum b value 

{bmin) ■ 



